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Abstract 

We examine the effective interaction of nonrelativistic fermions with an external vector field in 
superfiuid systems. In contrast to the complicated vertex equation, usually used in this case, we 
apply the approach which does not employ an explicit form of the pairing interaction. This allows 
to obtain a simple analytic expression for the vertex function only in terms of the order parameter 
and other macroscopic parameters of the system. 

We use this effective vertex to analyze the linear response function of the superfiuid medium 
at finite temperatures. At the time-like momentum transfer, the imaginary part of the response 
function is found to be proportional to Vp, i.e. the energy losses through vector currents are 
strongly suppressed. As an application, we calculate the neutrino energy losses through neutral 
weak currents caused by the pair recombination in the superfiuid neutron matter at temperatures 
lower than the critical one for the ^Sq pairing. This approach confirms a strong suppression of the 
neutrino energy losses as predicted in Ref. |2]. 
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I. INTRODUCTION 



Interactions among fermions in a superdense medium substantially modify their coupling 
to external fields. This problem is of a great importance both for a laboratory superconduc- 
tors and for a superfiuid baryon matter of neutron stars, where the thermal breaking and 
recombination of Cooper pairs is considered as a possible dominant mechanism of the neu- 
trino energy losses through neutral weak currents. The rate of reaction, in the ^5*0 superfiuid 
neutron matter, was for the first time estimated in Ref. Recent analysis [2i] has shown 
however that the approach used in this calculation contradicts the hypothesis of conserva- 
tion of the vector current in weak interactions and therefore considerably overestimates the 
neutrino energy losses. 

An alternative calculation of the same process was suggested in Ref. 3|], where the Green 
function technique was used to connect the rate of neutrino emission with the current-current 
correlation function, also called the retarded weak polarization tensor, and the corrections 
due to the strong interactions have been incorporated in the weak vertex function. The 
estimate was based on the idea that the dominant response of nonrelativistic baryons is due 
to the coupling of the baryon density to neutrinos. In this case the relevant input for the 
calculation is the imaginary part of the temporal component of the retarded vector-vector 
polarization tensor, Imll^^ {u, q) ,which was estimated by the authors in the kinematical 
domain 2A < u < oo and q = 0, to obtain Imll^^ {u > 2A, 0) ^ 0, where A is the superfiuid 
energy gap, and q = {u, q) is the transferred four momentum. It is easy to see, that this 
estimate contradicts the current continuity condition, uU^^ {u, q) = qiW^ {u, q), which tells 
us n°° {uj, 0) = when u > 2A. 

A one more calculation was recently done in Ref. {4] where the vertex renormalization is 
reduced to a study of the effective three-point vertices that sum-up particle-hole irreducible 
ladders in the scalar channel. The vector channel of the particle-hole interactions was 
parametrized by a constant (the lowest order Landau parameter) and the ideal summation 
is performed in the random phase approximation. The result of this calculation is also 
incorrect. Indeed, the whole longitudinal polarization function, as given by Eq. (35) and/or 
by Eq. (48) of this work, vanishes when the gap goes to zero. Thus, at temperatures higher 
the critical one for neutron pairing, both the real and imaginary part of the polarization 
function vanish instead of going over to those of normal (unpaired) Fermi liquid 
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A reasonable calculation of the medium response to external vector fields must satisfy 
the current continuity conditions. Neglect of this requirement leads to a considerable error 
in the estimate of the longitudinal response function. It is well known that the current 
conservation in superfluid fermionic systems is fulfilled in a nontrivial way. The interest 
to this problem has arisen many years ago in connection with the gauge invariance of the 
Bar deen- Cooper- S chr ieff er (BCS) theory 

B, a, a, a 

, [9|. But the problem is topical up 

to now fl, ml li- 
lt was realized that the current continuity equation can be satisfied if the interaction 
among quasiparticles is incorporated in the coupling vertex to the same degree of approxi- 
mation as the self-energy effect is included in the quasiparticle. This prescription yields the 
vertex equation explicitly incorporating two physical inputs - the pairing interaction and 
the energy gap which are not independent but connected by the g ap equation. Because the 
reahstic pairing interactions are momentum dependent [3], [l3], 14], particularly due to a 
short-range repulsive core, the above set of equations admits only numerical solution. 

On the other hand it is well known that the BCS theory of superfluidity needs the 
particular form of the pairing interaction only for a calculation of the order parameter (the 
energy gap A). The latter completely deflnes all the properties of a superconductor. Thus, it 
is reasonable to expect that, in the BCS approximation, the vertex function can be expressed 
only in terms of the order parameter, and other macroscopic parameters of the superfluid 
system, without of using of an explicit form of thepairing interaction. Up to now, the 
relevant solution of the vertex equation is obtained j?], only for the case when both the 
transferred energy and momentum are smaller than the superfluid energy gap, uj, |q| < A, 
and only for zero temperature, T = 0. In applications, however, we frequently deal with the 
time-like momentum transfer at flnite temperature T > 0. 

In this paper we focus on the analytic calculation of the corresponding effective vertex 
valid at flnite temperatures and for u and q, satisfying only u fi, |q| ^ pp, where 
H = p'jp/2M* and are the effective chemical potential and Fermi momentum of interacting 
fermions, respectively. Our calculation explicitly ensures the current conservation, using 
method that imposes this from the start. The essential idea of the method is as follows. 
We know that the external longitudinal fleld modifles the order parameter in the system. 
The corresponding self-consistent correction can be found directly from the condition of the 
current conservation. As was repeatedly discussed in the theory of superconductors, this 
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procedure is equivalent to the above summation of the vertex corrections [8j, [lOj. 

The vector vertex in the BCS approximation is to be considered as a starting point 
for incorporation of residual interactions in order to provide the current conservation in 
the superfiuid system. In turn, the residual interactions, evaluated in the random phase 
approximation are known to take into account the effects of the medium polarization (see 

n n 

e.g. [Ij], [15|), which renormalizes the sound velocity and the normal part of the vector 
vertex. This is beyond the scope of our consideration. 

The paper is organized as follows. Section 2 outlines some well known properties of 
normal nonrelativistic Fermi liquids at low temperatures and introduces the notations used 
in the following. In section 3, we derive the general expression for the effective vector 
vertex of electrically neutral fermions in the pair-correlated system and evaluate the obtained 
expression in the quasiparticle approximation. In section 4, we consider the linear response 
of the superfiuid medium in the vector channel. As an application, in section 5, we evaluate 
the neutrino energy losses through neutral weak currents caused by the pair recombination 
of neutrons in the crust of neutron stars. Section 6 contains a short summary of our findings 
and the conclusion. 

We use the system of units h = c = 1, and the Boltzmann constant fc^ = 1. 



II. NORMAL FERMI LIQUID 

In order to introduce the approach used in the following consideration this section out- 
lines some well known properties of normal nonrelativistic Fermi liquids at low temperatures. 
The degenerate fermion systems are commonly treated in the quasiparticle approximation. 
Within this framework, near the Fermi surface, the inverse Green function of normal (un- 
paired) fermions takes the quasiparticle form 

G^'ie,p)c^e-^p, (1) 

where 

2M* 2M* ~ M*^^ ^ ' 

is the energy of a quasiparticle relative to the Fermi level. The effective mass M* of a 
quasiparticle is connected to the Fermi momentum by means of the Fermi velocity Vp = 
Pf/M* which is small in a nonrelativistic fermion system. 
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In the lowest-order expansion ([T]) the wave-function renormahzation which accounts for 
the next-to-leading term in the expansion around the Fermi-energy is set to unity. For the 
same reason we neglect the imaginary part of the quasiparticle self-energy ~ \e\e. 

It is convenient to work in the particle-hole picture (known as the Nambu-Gor'kov for- 
malism) by introducing the quasiparticle fields as two-component (particle-hole) objects \I'p 
and using the Pauli matrices fi (i = 1, 2, 3) operating in the Nambu-Gor'kov space (see e.g. 
[3]). Then the Hamiltonian of the system of free quasi-particles takes the form 

Ho = Yl *J^P^3^. (3) 
p 

which corresponds to excited states in the particle-hole picture, while the ground state 
(vacuum) is the state where all negative energy "quasi-particles" (e < 0) are occupied and 
no positive energy particles exist. 

In terms of these fields, the vector vertex of a quasiparticle becomes the diagonal matrix: 

7> + ,.p)^f^""'^'''" ,° I. (4) 

\ -%{p + q,p) 

where the quasiparticle component is given by 

i,.{p + q,p) = ( 1,4- (p+^q 



^ ' M* V 2 

and the hole component is 7^ {p + q,p) = {—p, —p — q). 
Using the Pauli matrices this can be recast as 

7"=! = . (6) 

[liF(p+|q) if fi = 'i = 1,2,3 

The inverse quasiparticle propagator takes the following matrix form 

Go^ P) = ^ - ^p^3- (6) 



III. SUPERFLUIDITY EFFECTS 



In the particle-hole picture, the pairing represents a quasiparticle transition into a hole 
(and a correlated pair). Therefore the self-energy arising due to the pairing interaction has 
the off-diagonal form in the Nambu-Gor'kov space. In the absence of external fields this 
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anomalous self-energy is real-valued and can be written in terms of the energy gap A arising 
in the quasiparticle spectrum. Then the BCS inverse propagator of a quasiparticle has the 
following form (see e.g. (9|): 

G-i(£,p)=£-epf3-Afi, (7) 
The the effective vertex in the pair-correlated system can be written as 

where 7^ is the "normal" vertex, as given by Eq. (jl]), and S*^^-* is the linear correction to 
the anomalous self-energy of a quasiparticle in the external vector field V^. Apparently this 
correction is due to varying of the energy gap in the external field. By taking the new gap 
in the form 

A = A + A^^) {V^) 

we can write 

y A*(i) (v^) 

Now we use the fact that, in the case q -C p^, the linear effect of the weak field pertur- 



bation is to change the phase $ (V^) but not the magnitude of the gap parameter j8|, 10|. 
The field- induced correction to the amplitude of the gap is proportional to A/yU ^ 1 and 
thus can be neglected. We therefore write the energy gap as A = Aexpi$ (V^), where A is 
the (real) gap in the absence of external field. Apparently the linear correction to the gap 
A = A — A is proportional to the external field and can be written as 

A(i) (V) = A • (e'*(^^) - 1) ~ 2A • $ (\/^) ~ -Q^ (q) (q) , 

^y^.^ _ ^ . (^^-^'^>iv-) _ 1) ~ • <!> {V) ^ (g) (g) , (9) 

where (g) is the unknown kernel connecting the correction to the energy gap with the 
external weak field. Then the effective vertex becomes 

f> + ,;rt=(^'"'^'''" -'^^^'^ ], (10) 

and the problem reduces to calculation of the unknown vector function (tu, q) . 
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The relation between the modified vertex F'^ and the quasiparticle propagator G is given 
by the matrix Ward identity 

q^-V, {p + q,p) = {p + q)h- %G-' (p) . (11) 

Applying this identity to expressions ffTOj) and ([7]) we obtain the following well known relation 

QuQ" (g) = -2A. (12) 

For further progress, let us consider the medium response Il'^'^{q) to the external vector 
field. This calculation needs the quasiparticle propagator which, in contrast to the inverse 
propagator ([7]), must be uniquely defined at finite temperature. Since the matter is assumed 
in thermal equilibrium at some temperature, we employ the Matsubara calculation tech- 
nique. In this case the quasiparticle propagator can be obtained by inverting Eq. (I7l), and 
substituting e = in {2n + 1) T with n = 0, ±1, ±2.... Then the retarded polarization tensor 
1^*^(5) is given by the analytical continuation of the following Matsubara sums (see e.g. 

U^'^iUn,, q) = T22 J (^^^'^^ [^^ +U}m,P + q.;Pn, P) G {pn + U}m,P + q) 

Pn 

X f'' {pn + u;^, p + q; p) G (p„, p) (13) 

onto the upper-half plane of the complex variable u. Here Pn = tt {2n + 1) T, and 
Um = 27rmT, with m, n = 0, ±1, ±2..., are the fermionic and bosonic Matsubara frequency, 
respectively; and the quasiparticle propagator has the following form 

GIp.pjJ^'^-P) ^•(""•P' V (14) 

\F{p^,p) -Gt(p„,,p) J 

The particle and hole components are given by the on-diagonal elements, with (p„, p) = 
G {—pn, — p) while the off-diagonal elements, (p„, p) = F {—pn, — p), represent the anoma- 
lous contribution caused by the pairing interaction. According to Eg. ([7]), in the BCS 



approximation the components are given by the following expressions 17l |: 



— ^Pn — ^ A 

pI + pI + £p 

Ft (p„, p) = F (p„, p) , Gt (j9„, p) = (15) 
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where is the energy of a quasiparticle 

= ^/^^^■ (16) 

It is important to notice that the magnitude of the energy gap in the above expressions 
depends on the temperature, A = A (T). We assume that this function is found from the 
corresponding gap equation. 

Using expressions (fT^ . (|TOl) in Eq. (|T3l) gives 

ii,Aq) = KAi) + ^,{q)Qu{q). (17) 

where the one-loop integrals (c<j, q) and (c<j, q) are defined as the analytical continua- 
tion of the following Matsubara sums: 

^^u (t^m, q) = T / ——3 [itiluG {pn + UJm, P'+Q) G p') 

+ llllG^ {Pn + ^m, p'+q) G^ {p 

- bf^ll + lllu) F (p„, p') F {pn + ^rn, p'+q)] (18) 
K {^rn, q) = T / ——3 b/xG (Pn + ^m, p'+q) F {pn, p) 

Pn ^ ^ 

+ llF {pn + LUm, p'+q) G^ {prr, p') 
- llG^ (Pn + UJm. p'+q) F {p 

n; P j 

-1^,F {pn + uJm, p'+q) G {pn, p')] • (19) 

The unknown vector function Q^, {q) can be found from the requirement that the polar- 
ization tensor (IT7|) obeys the current conservation, Ti^'^q^ = 0, and qu^^^ = 0. The first of 
these relations is satisfied automatically if the effective vertex T ^ obeys the Ward identity 
f fTTj) . the second one is to be satisfied by a proper choice of the function (g). 

With the aid of relation (fT2|) the above conditions can be written as the two coupled 
equations 

A'^^(g)g,-2A A^(g) =0, (20) 

q,k-^{q) + qA''{q)Q^{q) = Q. (21) 

Using the symmetry property, q^K^^ = Iv^^qu, we find 

2A 



Inserting this in Eq. (fTOj) we obtain the effective vector vertex in the following form: 

Thus the BCS approximation reduces to the calculation of the vector function A'^ (g). In 
Eqs. ( fTSj) . ( fT9|) . the summation over the Matsubara fermionic frequency can be performed 
exactly. The subsequent analytical continuation yields 



Aoo (t^, q) 



2 (27r)^ 



_ ( 1 + ^p^p+q 1 $_ 

^p^p+q 



(24) 



Ao (-, q) = - / r^T^ + , (25) 
A. q) = / f P.+iq.l ($, + $-) . (26) 



M* J 2 (27r)' V 2 



^p^p+q 



\2 



,.A.(.,q) = -/-^^%i^ + (27) 

where the following notations are introduced: 

^ £p+q±£p / ^ £pj^\ _ 

(ep+q ± -{UJ + zOf \ 2T 2T J ^ ^ 

Eq. (!23|) is the general result valid for arbitrary momentum transfers. By the direct sub- 
stitution one can verify that this effective vertex satisfies the Ward identity (|TT|) . providing 
the current conservation in the system. 

As we can see, in the BCS approximation, only the longitudinal components of the 
vertex are modified because the transverse components of Aj vanish at the angle integration. 
Making use of this fact we need to calculate only the temporal component of the vertex 
because the remaining components can be found from the current continuity condition. 

General expression for the vertex function, as given by Eqs. (123|) - (128|) . has a complicated 
form. Significant simplification is possible, however, due to the fact that the quasiparticles 
are nonrelativistic, i.e Vp <^ 1 and because we are interested in q ^ pj?. The latter condition 
implies q/M* = Vpq/pp <^ Vp and we may neglect the quasiparticle recoil. Indeed, in 
the cases typical for astrophysical applications one has Vp ~ 0.1 while q/M* ~ T/M* ~ 
10"'^ ^ 10^'^. Under these conditions the contributions, caused by the recoil are estimated 
as qV {M*u) ~ (q/M*) (q/tu) ~ Vp (q/p^) (q/iu). 



We have 

- = q Vf cos ^ + ^ = q Vp (^cos 6 + (29) 
where 6 is the angle between p and q momenta. Insertion in Eq. (12 7p gives 

g.A.(u;,q) = -qVJ / Lsg + -^y fP f ($+ + $-)■ 
J \ ^Pf/ (27r) ^£^p£^p+q 

In this expression, contributions due to the odd powers of cosine vanish after angle integra- 
tions. To the lowest order in small parameters we may take the integrand at q = to find 
the following relation 

1 V'^ 

QiAi {uj,q) ~ -q'^Ao (w,0) 

UJ 6 

valid up to accuracy Vp. Apparently to the same accuracy we can write 

QiAi {u, q) ~ -q^c^Ao {uj, q) , (30) 

UJ 

where = Vp/y/S is the sound velocity in the Fermi gas. 
Insertion in Eq. (123!) immediately gives 

f o(p + q;p) ^ ra - 2^Af2 , , , , (31) 
The space component of the longitudinal vertex can be found from the Ward identity: 

2 2 

qf ^ ^p+q - - 2zAf2 (32) 

It is necessary to stress that the above expressions are obtained from the exact expression 
Eq. (!23|) to accuracy Vp <^ 1 valid in the nonrelativistic system. To this accuracy the 
temperature dependence of the effective vector vertex is restricted to the gap function A (T). 

It is interesting to notice that the approximate Eqs. fl3ip and fl32l) reproduce formally 
the vertex functions as obtained in Ref. j3] for the case of small transferred energy and 
momentum, u;, q ^ A, and zero temperature. In contrast, our analysis shows that Eqs. 
( |3T|) . (!32|) are valid under considerably weaker conditions uj <^ jj, = p|,/2M*, q <^ p^? and for 
a time-like momentum transfer, |q| < cj, as well. The energy gap A = A (T), in Eqs. fl3Tl) . 
and (1521) . depends on the temperature below the critical temperature T^. At temperatures 
T > Tc the gap vanishes. In this case, Eqs. fl^ . and fl5^ give the vector vertex (jS]), as 
obtained for normal Fermi liquids. 

As is well known, the poles in the vertex function at cu^ = q^c^ indicate the existence 
of the collective acoustic-like excitations in the condensate propagating with a velocity Cg. 
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The second term in (13T|) . and (132|) is the result of the couphng of % to the collective mode. 
This can be understood in the following way. f ^ contains matrix elements for creation or 
annihilation of a pair out of the vacuum. This process can go through the virtual collective 
intermediate state. This collective contribution plays the important role in the conservation 
of the vector current in superfiuid systems. 



IV. LINEAR RESPONSE IN THE VECTOR CHANNEL 



Having determined the effective vertices, we turn to evaluation of the complete polariza- 
tion tensor, which according to Eqs. f|T7j) . fl22|) is given by 

2A 



n^^ (g) = Ami. (g) - 



(g) K (q) , 



(33) 



^aA^ (g) 

With the aid of relation fl20l) this expression can be transformed to the form which explicitly 
exhibits the continuity of the current: 

1 



(g) = A^^^ (g) 



gVA^„ (g)A^, (g) . 



q^q^Axs (g) 

Indeed, this form obeys li^yq^ = q'^H^u = 0. Using this fact we decompose the polarization 
tensor (|33|) into the sum of longitudinal (with respect to q) and transverse components 

n''^ (g) = (g) (^1, ^n^ ' (^1, ^n^ + (g) g^' {S'^ - n^n^) g^\ (34) 

Here n = q/q is a unit vector; the longitudinal and transverse polarization functions are 
defined as 

H, (g) = H°° (g) , H, (g) = ^ {6'^ - n^n^) ff^' (g) . (35) 

As it follows from the definition (fT8|) . A^j, (cj, q) represents the polarization tensor in 
the so-called one-loop approximation. The transverse external field does not influence the 

171], therefore the transverse polarization function 



anomalous self-energy of a quasiparticle 



is given by the one- loop expression (!36|) . The calculation yields the well known expression 



(^,q) 



1 



+ 1 1 



d^P 2 ■ 2n 

o p sm 

{2'nf 



+q 



PS,p+q 



(36) 



Compare this to the linear response of superconducting electrons to external transverse 



electromagnetic field 



17| 
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In contrast, the longitudinal external field varies the order parameter in the system. This 
physically manifests itself as the collective acoustic excitation in the condensate. Therefore 
the longitudinal polarization function has the additional contribution, which takes into ac- 
count the collective motion of the condensate providing a continuity of the current in the 
system. We obtain 



2A 



Ao (g) Ao (g) 



^aA^ (g) 

According to Eq. (l30|l . up to accuracy Vp this can be simplified as 

2Aw 



Ui (^,q) ~ Aoo (w,q) 



or, identically. 



+ 



2(27r) 
up' — q^c^ 



1 - 



PSp+q 



2(27r) 

By the use of the following trick 



a;2 — q2c2 



Ao (^,q) 



1 + 



p'^p+q 



^p^p+q 



uo"^ — q^c^ 



1 + 



2 - q2c2 



we can write 



dh 



2(2 



+ 



q^c2 



C<j2 — q2c2 



1 - 



PS,p+q 



A^ 



^p^p+q 



2A2 



$-1 



^p^p+q + 



^p^p+q 



(37) 



(38) 



2 (27r)^ £^p£^p+q 

and substitute ^p+q as given by Eq. ( l29i) . Then with the aid of the series expansion of the 
integrand up to accuracy (qVp/A)^, and the averaging over angles after some simplifications 
we arrive at the following result: 



+ 



X 



2T2 



q^c2 



rfp(l-l^ln " + ^"^ 
2 qup uj — qup 

A2 



Pi 



[uj + ^0)' - q2c2 27r2 
q2c2 1 



dp— tanh — 



(cu + zO) -q2c2 2T27r2 
1 u; UJ + qUr 



dp- 



A2 



2qup 



In 



UJ 



tanh^ ^) 
2TJ 



tanh^ ^ 
2T 



(39) 
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where 



Up = Vf—- 



p 



is the velocity of a quasiparticle. 

As given by this expression, contributions into the longitudinal polarization of a superfluid 
fermion system arise both due to the motion of quasiparticles (the first line) and the collective 
motion of the condensate (the second line). There is also the mixed term given in the last 
two lines of Eq. (l39l) . 

At T > Tc the energy gap vanishes and Eq. fl39l) reduces to the well known expression 
for the longitudinal polarization of a normal Fermi gas: 



The quasiparticle contribution and the mixed term vanish at zero temperature, when all 
the fermions are paired and no quasiparticle excitations exist. Therefore, at zero tempera- 
ture, the longitudinal polarization arises only due to the collective motion of the condensate. 
We obtain 



A. Imaginary part 

The imaginary part of the polarization function deserves of a special consideration because 
it contains a complete information on the energy losses and scattering of external particles 
in the medium. Here we have to distinguish the space-like, > cu^, and time-like, < u"^, 
kinematical domains. One can easily see that the polarization function ( l39i) has no imaginary 
part in the time-like kinematical domain. This due to the fact that the polarization function, 
as given by Eq. (!39|) . is calculated to accuracy Vp. As we shall see, in the time-hke 
kinematical domain, the imaginary part is proportional to Vp. 

1. Space-like momentum transfer 

In the space-like kinematical domain, q^ > cj^, the longitudinal polarization function has 
the imaginary part due to the pole contributions of the collective mode at u"^ = q^c^ and due 
to the Landau damping at u; < qup. The corresponding expressions can be readily found 




n/ (g;T = 0) 



PfM* 




(40) 
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from Eq. (139|) . At finite temperature this expression has a cumbersome form and represents 
no practical interest in the nonrelativistic case. Therefore we shall restrict our consideration 
of the space-like momentum transfer to the case of zero temperature. From Eq. (l4Up we 
readily find 

Imn, (g; T = 0) = -^q^c^5 {u' - q^c^) signer 

Since no quasiparticles exist at zero temperature, the imaginary part arises only due to 
excitations of the collective motion of the condensate. This form satisfies the /-sum rule 

3 

2 



f ujlmUi(q;T = 0)duj = 



2. Time-like momentum transfer 

We focus now on the imaginary part of the retarded polarization function at the time- 
like momentum transfer, < cj^, important for applications. As mentioned above, this 
imaginary part needs the calculation of corrections of the order Vp. 

At > 2 A and q < uj the imaginary part in Eqs. fl24l) - fl27j) and fl36|) arises from the pole 
of the function $_|_ uo = ^p+q + ^p. Calculation of the imaginary part of the transverse 
polarization function fl36|) is trivial. We get 



Im n, {u, q) = — r>^M*^!^^|i£^ tanh — + o( V^^, Vl) (41) 

where B (x) is the Heaviside step function. 

Calculation of the imaginary part of the longitudinal polarization function is more com- 
plicated. First of all, we evaluate the one- loop integrals up to accuracy Vp. This gives: 

Aoo(cc', q) ^ ao{uj, q) + a^iuj, q) V| + O {v^^^, V^^ , (42) 

Ao(^, q) - boiiu, q) + hiiu, q)l^| + O (v^^, V^^ , (43) 
qM^, q) - C2(^, q)^| + O (v^^, , (44) 
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where 



PfM* ^ uj 

Iman = =tanh— — , 

vr ujVu^ - 4A2 4T' 

Imao = ;t = tanh — — 

1 - tanh^ —uj]], (45) 



cj2-4A2 16T2 V 4T 



Refen = n — /" de , — rtanh—— , 

- A2 4^2 - 4T' 

T A PfM* a ^ 

im On = = tann — -, 

2n Vcu2 - 4A2 4T' 

PfM* q2 A ^ 

imoo = , tann — - 

271 3a;2 Vu;2 _ 4^2 4T 

/2A2 + CU2 (c^2_4^2) _ ^ ^ ^ ^ 



Rec2 = — q^ ^ , V / & , :— rtanh — 

3^ 7r2 ^^2 _ ^2 4^2 _ ^2 42^ 



Imc2 = — -^^^ — =tanh— — ■ (47) 

TT 3uj Vcj2 - 4A2 4T ^ ^ 

In the above, symbol V means the principal value of the integral. 

As already mentioned, at q2 < ti;2^ ^]-^g imaginary part of 11/ is small as compared to the 
contributions evaluated in previous sections to accuracy Vp. Therefore, instead of Eq. (!38ll 
we start from the exact form of the longitudinal polarization function, as given by Eq. ([3j 
in order to obtain the following series expansion over powers of Vp 

A \ / A /, 1 \ \ _o Ad 



ao - 2-60 ) + ( a2 - 2- ( 62 + -Co ) ) - 2^?^^- 

UJ / \ UJ \ UJ UJ'^ On 



Thus we have 



Im Hi ~ ( Im ao — 2 — Im bo 

UJ 



+ (lma2 - 2-lmh - 2^lmco] - V^^Imf^C2^ ... 
V ^ uj^ J uj-^ \bo J 

Upon substituting Eqs. fHB]) . (H7|) in the first and second terms of this expression we 

find that these leading terms mutually cancel 



A 

Imao — 2 — Im6o = 0, 

UJ 



15 



A A 
Im a2 — 2 — Im 62 — 2— Im cq = 0, 

in accordance with the results obtained in previous sections. 
In the third term we find identically 

^ ^ 2g^ 

and thus 

A g2 
Imllj ~ -4—— — rlmc2 

We finally obtain the following lowest-order expression 

I.nn,(..q).lc;p,Ar^^ff';-^f tanhi^ (48) 

TT cu^v^^ — 4A^ 4i 

valid at q < UJ. 

We see that, in comparison with the one-loop approximation, as given by Eqs. P2|) 
and (l45l) . the self-consistent result (148|) is suppressed by many orders of magnitude. The 
suppression factor is ~ 10~^. 

V. NEUTRINO EMISSION AT THE PAIR RECOMBINATION 

As an application of the obtained results we consider the neutrino-pair emission through 
neutral weak currents occurring at the recombination of quasiparticles into the condensate. 
The process is kinematically allowed due to the existence of a superfiuid energy gap A, 
which admits the quasiparticle transitions with time- like momentum transfer q = (ti;,q), as 
required by the final neutrino pair. 

We consider the total energy which is emitted into neutrino pairs per unit volume and 



time which is given by the following formula (see details e.g. in 15l|): 

where Gp is the Fermi coupling constant, is the neutrino weak current, and ^'^^^^y- is 
the retarded weak polarization tensor of the medium. The integration goes over the phase 
volume of neutrinos and antineutrinos of total energy uj = uji + 002 and total momentum 
q = q^^ + q2. The symbol Ylu indicates that summation over the three neutrino types has 
to be performed. 
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We shall consider the ^5*0 pairing of neutrons, which takes place in the superfluid crust 
of neutron stars. In this case, the total spin of a bound pair is zero and the axial-vector 
weak interaction occurs only due to small relativistic effects which are usually omitted in 
the calculations [1]. Therefore we concentrate on the vector part of the weak interactions. 
In the vector channel, the weak polarization tensor is given by n^^g^k i^) ~ CyTl^'^ (g) , where 
Cv is the vector weak coupling constant of a neutron, and 11'^'^ (g) stands for the polarization 
tensor as given by Eq. flM|) . 

By inserting j d'^qS'^^^ il ~ Qi ^ Q2) = 1 in this equation, and making use of the Lenard's 
integral 

/ mi^''^ - - = Y (^^'^^ ~ ® ® ' 

where 9(x) is the Heaviside step function, and g^i, = diag(l, — 1, — 1, — 1) is the signature 
tensor, we can write 

487r* Jo Jo exp [^) - 1 

where A/^, = 3 is the number of neutrino flavors. 

Due to the current conservation g^g^^Imll^'^ (g) = 0, an we flnd 

Q = l^AT. r du r dq ^"'~f " [(^ - 1) Imn. (g) + 2Imn, (g) 
487r4 Jo Jo exp (f) - 1 LVq / 

where the functions Imllj j (g) are given by Eqs. (HTl) and (HHI) . 

After simple calculations the neutrino energy losses are found to be: 

-^.x (50) 



where y = A/T, z = \/ x'^ + y'^. 

We see that the neutrino radiation via the vector weak currents in the nonrelativistic 
system is suppressed by several orders of magnitude with respect to that predicted in 
The suppression factor is ~ V^. This result is consistent with the result obtained in Ref. 
[2], by the use of the Fermi Golden rule. 

VI. SUMMARY AND CONCLUSION 

We have developed a unifled approach for calculating the longitudinal response function of 
superfluid fermion system at flnite temperatures. The effective vertex in the vector channel 
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is analytically expressed in terms of the order parameter, and other macroscopic parameters 
of the superfiuid system, without of using of an explicit form of the pairing interaction. 
The general expression for the vertex function is given by Eq. (1231) . When evaluated to 
accuracy Vp <^1, this expression becomes very simple and takes the form, as given by Eqs. 
( 1ST]) . ([22D- The latter formally reproduce the result obtained in Ref. [7j for the case of zero 
temperature and small transferred energy and momentum. In contrast, our result is valid at 
finite temperatures and for u and q, satisfying only u <^ /^,|q| <^ Pf- As we found to this 
accuracy <^ 1 the temperature dependence of the effective vector vertex is restricted to 
the gap function A (T). This justifies the approach used in Ref. where this form of the 
effective vertex was used in the calculation of the neutrino energy losses. 

Special attention was paid to preserving the Ward identity for the vertex function and 
the relations for the polarization tensor components, that are implied by the conservation of 
the current. The transverse and longitudinal polarization functions, as represented by Eqs. 
( !36l) . ( !39l) . are calculated up to accuracy Vp. 

The imaginary part of the response functions, as given by ( l36l) . ( l39l) . arises due to the 
pole contribution of the acoustic-like collective mode and due to the Landau damping. Both 
contributions vanish at the time-like momentum transfer, uj > 2A and q < uj. In this 
kinematical domain, the imaginary part of the polarization functions is due to breaking of 
the correlated pairs and recombination of the quasiparticles back into the condensate. The 
corresponding analytical expressions are found to be proportional to Vp, as given by Eqs. 

As an application of obtained results we have calculated the neutrino-pair emission 
through neutral weak currents occurring at the recombination of quasiparticles into the 
condensate. We found the rate of neutrino energy losses in the vector channel is suppressed 
as compared to the one- loop results by a factor Vp. The magnitude of the suppression is in 
accordance with the one predicted in Ref. . Thus the neutrino emission from the singlet- 
correlated neutron matter is mainly due to the axial- vector currents. As is well known 
the corresponding neutrino energy losses are proportional to Vp and thus also small. Ap- 
parently, the modifications to the neutrino emission rate through the pair recombination 
process, as found above, call for a detail reassessment of their role in the late-time cooling 
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of neutron stars. 
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